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^ l' Abstract. We study two-flavor collective neutrino oscillations in the dense-neutrino region 

D , above the neutrino sphere in a supernova (SN). The angular dependence of the neutrino- 

neutrino interaction potential causes "multi-angle" effects that can lead either to complete 
kinematical decoherence in flavor space or only to small differences between different trajectories. 
This nonlinear system switches abruptly between "self-maintained coherence" and "self-induced 
^ ' decoherence" among the angular modes, depending on the strength of the deleptonization flux. 

I For a realistic SN the quasi single-angle behavior is probably typical, simplifying the numerical 

■ treatment and probably allowing for the survival of observational features of flavor oscillations. 

(N 

The neutrinos streaming off a collapsed supernova (SN) core are so dense near the neutrino 
! sphere that they produce a significant refractive effect for each other, leading to collective 

' oscillation effects. The practical importance of these nonlinear phenomena was only recently 

recognized and studied in a series of papers [IlElEllliaiSElElElIIOlIIIlIiaiiailllllS]. An 
introduction was given at this conference by Eligio Lisi |16j . 

■ The SN neutrino fluxes are thought to obey the hierarchy Fi,^ > Fp^ > Fj^^ where stands 
5^ , for any of z^^, Vt, and Vr- In other words, there is an excess of feZ^e pairs relative to Vx^x- In 

an inverted hierarchy situation, this pair excess converts collectively into Vx^x pairs, a process 
that does not violate flavor-lepton number and thus does not require mixing: it could also 
proceed as an ordinary pair annihilation process. Neutrino refraction causes this pair process to 
proceed collectively and very fast, almost independently of the mixing angle, i.e., we have to do 
with a collective "speed-up effect" [U [2]. The unpaired Ve excess flux from deleptonization 
is conserved [5]. In the adiabatic limit, it is the low-energy part of the spectrum that 
survives, leading to a step-like feature in the Ve spectrum (a "spectral split" [TO] caused 
by a "step-wise spectral swapping" [U dH). All of these effects happen in the dense- neutrino 
region that typically extends from the neutrino sphere out to a few hundred kilometers. If the 
ordinary MSW resonances occur at larger radii, the collective phenomena and subsequent MSW 
transformations are independent, the former producing the initial condition for the latter. For an 
inverted hierarchy the only effect of matter in the collective-transformation region is a decrease 
of the effective mixing angle [5l[8]. Conversely, if the matter density profile is very shallow, the 



X 



ordinary MSW effect can prepare tfie initial condition for the collective effects |14j . 

Mixed neutrinos are described by matrices of density pp and pp for each (anti)neutrino mode. 
The diagonal entries are the usual occupation numbers whereas the off-diagonal terms encode 
phase information. The equations of motion are idtQp = [Hp, Qp], where the Hamiltonian is [TT] 

Hp = 17p + V + V2 ^ (£>q - ^q) (1 - Vq • Vp), (1) 

Vp being the velocity. In the mass basis, the matrix of vacuum oscillation frequencies is 
Op = diag{m1,m2,m'^)/2E. The matter effect is represented, in the weak interaction basis, by 
V = GpUB diag{Yf,, 0, Y^^). For antineutrinos the only difference is Op —Op. The effective 
tau-lepton density Y^^ ~ 10~^ arises from radiative corrections [18] and can be important in a 
genuine three- flavor treatment of ordinary [19] or collective [15] SN neutrino oscillations. For 
the latter, the influence of Y^^ can be rather sensitive to deviations from maximal 23-mixing. 

The angular factor (1 — Vq-Vp) in Eq. ([T]) derives from the current-current nature of the weak- 
interaction Hamiltonian. In an isotropic ensemble it averages to unity, whereas the neutrinos 
streaming off a SN core are strongly non-isotropic so that different angular modes experience a 
different strength of the neutrino-neutrino interaction potential. One may expect that this effect 
leads to kinematical decoherence among angular modes and thus to flavor equilibrium [HEj. In 
a symmetric ensemble with equal densities of neutrinos and antineutrinos, this is indeed the 
case. Such a system is highly unstable in that an infinitesimal deviation from exact isotropy is 
enough to trigger an exponential run-away towards flavor equilibrium On the other hand, a 
numerical simulation of the flavor evolution of SN neutrinos revealed that multi-angle effects were 
small: All angular modes evolved nearly collectively, very similar to an isotropic ensemble [1]. 
Likewise, in a multi-energy system every energy mode feels a different Hamiltonian, yet the 
evolution is collective, i.e., multi-energy effects do not lead to kinematical decoherence. 

We have performed a numerical exploration of multi-angle effects in a spherically symmetric 
system where the neutrinos are emitted from a "neutrino sphere" [8]. We have identified the 
asymmetry between the and flux as the crucial parameter. Realistic deleptonization fluxes 
in SNe seem sufficient to suppress multi-angle decoherence. Therefore, in practice multi-angle 
effects seem to be a subdominant feature of collective SN transformations. 

A better understanding of collective oscillations can be developed in the two-flavor case in 
terms of the usual flavor polarization vectors Pp to express the matrices pp. The equations of 
motion are (?fPp = Hp x Pp. Ignoring the ordinary matter effect, the "Hamiltonian" is 

Up=U;B + V2GfJ ^ (Pq - Pq) (1 - Vq • Vp), (2) 

where B is a unit vector in flavor space in the "mass direction" and to = Am'^/2E the vacuum 
oscillation frequency. In an isotropic situation where the velocity terms average to zero, the 
second term is of the form p, D where the vector D is the difference between the total neutrino and 
antineutrino polarization vectors and p = ^/2Gp Uy the neutrino-neutrino interaction strength. 
The individual Hp = wB -|- all lie in a single plane. In the adiabatic limit the polarization 
vectors follow the Hamiltonians so that they also lie in this co-rotating plane. This observation 
explains the collective nature of the multi-energy evolution (all polarization vectors stay in the 
same plane that rotates around B with a certain frequency Wc) and explains the spectral splits 
in that the final Hp in the co-rotating plane are either aligned or anti-aligned with B [9l [10] . 
This picture also illustrates a crucial difference to the symmetric system where initially D = 
and the overall evolution is purely pendular [5] — there is no co-rotating plane. 

In a spherically symmetric situation every mode is characterized by its vacuum oscillation 
frequency uj and its angle 9 relative to the radial direction, providing v = cos 6 as the radial 



velocity. Now the individual Hamiltonians are Hp = ti;B+/i(D— t;pF) where F is the flux term of 
the neutrino ensemble, i.e., the same as D, but every mode weighted with Vp. If all polarization 
vectors lie in a single plane, then also all Hp are in that plane so that an evolution in a co- 
rotating plane is self-consistently possible. Our numerical simulations show that those cases 
with little kinematical decoherence correspond to the polarization vectors essentially staying 
in a co-rotating plane with some zenith- angle spread. On the other hand, the vectors B, 
D and F do not have to stay in a single plane and in fact, this appears to be an unstable 
arrangement. Kinematical multi-angle decoherence corresponds to strong deviations from this 
coplanar situation. In a spherically symmetric situation, the neutrino-neutrino interaction 
strength decreases with r~^. In a toy model with an artificially slow decrease of /i(r) one 
can make the evolution arbitrarily slow and adiabatic. In this case multi-angle decoherence 
appears to be unavoidable. 

It appears that in a realistic SN the collective evolution is slow enough to be essentially 
adiabatic with regard to the development of a spectral split, but fast enough that the co-planar 
arrangement of the polarization vectors survives. This picture may provide the key to an analytic 
understanding of the conditions for multi-angle kinematical decoherence. 

The assumption of spherical symmetry of the overall system severely restricts possible 
solutions. The evolution is a one-dimensional problem along the radial direction. It remains 
to be studied if this symmetry has an important impact on kinematical decoherence, i.e., if a 
system with fewer symmetries would decohere more easily. 
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